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Recent progress in apphcation of higher order QCD calculations to jet and 
inclusive particle production in photon induced collisions is reviewed. At- 
tention is paid to theoretical uncertainties of such calculations, particularly 



I those coming from the choice of renormalization and factorization scales. 



PACS numbers: ll.lO.Hi, 12.38.-t, 12.38.Bx, 13.85.Ni 

1. Introduction 



From the variety of hard processes involving initial state photon I will dis- 
cuss jet and inclusive particle production in 7p and 77 collisions with em- 
phasis on recent phenomenological applications. I will concentrate on the 
discussion of theoretical uncertainties because they are often bigger than 
experimental errors thereby complicating the interpretation of data. 
(— I I Let me start with brief recollection of basic facts concerning the structure 

and hard collisions of protons and photons. For details see Parton 
distribution functions (PDF) of the photon depend on the factorization 
^ ' scale, denoted M below, and satisfy the system of evolution equations 



dS(x,M) 
dlnM2 

dG{x,M) 
dlnM2 

dlnM2 



Sskq + Pgg ® S + PgG ® G, (1) 
kG + PGq^^ + Pgg ® G, (2) 
5mkq + Pns ® qm, (3) 



where M), q^s and G{x, M) stand for standard singlet, nonsinglet and 
gluon distribution functions and the splitting functions Pij and ki are given 
as power expansions in as{M). For the proton the inhomogeneous splitting 
functions kg and kc are absent. 

General solution of the evolution equations can be written as the 
sum of a particular solution of the full inhomogeneous equations and the 
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general solution of the corresponding homogeneous ones. For instance, the 



solution of © including only the lowest order terms k^^ and Pqq^ and 
vanishing at some factorization scale Mq, is given in momentum space as 



,(0) 
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a-Nsin), (4) 
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2vr/3o 1 - 2P,V°\n)//3o 



(5) 



It is often claimed that because of the presence of ag in the denominator of 
(dJ) quark distribution functions of the photon behave as a/us- As argued 
in [2] this is misleading and PDF of the photon actually start at the order a. 
The standard assignment of the order of PDF of the photon has important 
implications for the definition of finite order approximations, but in this 
talk I will not pursue this point further. 

The cross section of inclusive hadron h production in 7p or 77 collisions 
has the generic form 

a{-fB -^hX)ocY^ Da/^i{M^) ® D^,/b[Mb) ® a{ab ^ cX) ® D^{Mh) (6) 

a,b,c 



where the sum runs over quarks, gluons and photons (in the direct photon 
channels) in the photon or proton. The hard partonic cross sections can be 
expanded in powers of as as follows 



aiab ^ cX) = a1{ii)a'''' {ab ^ cX) + a^+^(/i)c7''^^(a6 cX) + ■ ■ ■ (7) 

with k = 1 for direct photon processes and k = 2 for resolved ones. For jet 
production fragmentation functions are replaced by some jet algorithm. 



2. Remarks on phenomenology 

In (prri) I have distinguished the factorization scales of the beam particles 
from the fragmentation scale Mf^ and the renormalization scale //. The 
latter appears only in perturbative expansion of the partonic hard scatter- 
ing cross sections ((TJ. There is no theoretical reason to identify any two 
or more of these four scales as they come from the treatment of different 
singularities. Unfortunately, in most phenomenological analyses all the four 
scales are set equal and identified with some "natural physical scale" Q: 
fi = My = Mp = Mfi = M = Q. However, choosing the renormalization 
scale does not fix Os because Osin) depends beside the renormalization scale 
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also on the renormalization scheme (RS). Consequently, the same choice of 
the renormalization scale gives different results in different RS! The choice 
of the RS is in fact as important as that of renormalization scale /i, but 
there is no "natural" RS! The conventional procedure of working in MS RS 
is entirely ad hoc. Choosing the scales and schemes should be based on 
the investigation of the dependence of finite order approximants on renor- 
malization and factorization scales as independent parameters and should 
reflect the possible existence of regions of local stability. 

In comparing data with theoretical calculations the latter are quoted 
with estimates of their "theoretical uncertainty". This quantity usually 
comprises several components: dependence on PDF and FF, hadronization 
effects and scale and scheme variation. Hadronization corrections are usu- 
ally claimed to be small, but as we shall see in the case of jet and inclusive 
particle production in 77 collisions it does not have to be always the case 
and further investigation of this point is certainly worth the efforts. 

Whereas the uncertainty reflecting the dependence on PDF, FF and 
hadronization is reasonable well-defined, the conventional way of estimating 
the uncertainty due to the freedom of choice of scales and scheme is definitely 
not. Identifying all scales with some "natural scale" Q and varying this 
common scale M typically in the interval Q /2 < M < 2Q makes little sense. 
First, the results still depend on the selected renormalization scheme, and, 
second, the chosen range of multiplicative factor (1/2,2) is again entirely 
arbitrary. 

Another important aspect of the comparisons of theory with data con- 
cerns the relation between LO MC event generators, like HERWIG, RAP- 
GAP and PYTHIA, and NLO parton level calculations. Some of the features 
of full NLO QCD effects are mimicked within LO MC event generators by 
means of parton showers and noncollinear kinematics of initial state parton 
emissions. Moreover, LO MC use different input (PDF, FF and Og) which 
were extracted in LO global analysis from data, and so have a chance to 
describe also other data. 

3. Jet production in 7p and 7*p collisions 

Theory and phenomenology of jet production in (quasi)real 7p collisions 
has recently been reviewed in [3]. Extensive application to HERA data has 
lead to good agreement, except for some IR sensitive quantities but suffers 
from non-negligible scale dependence of existing QCD calculations. 

In this talk I will concentrate on the phenomenological application of 
QCD calculations to HERA data in the kinematic region of moderate Q^, 
where the transition region between photoproduction and deep inelastic 
scattering takes place. This region is of particular interest because of the 
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Fig. 1. Comparison of HI data on forward jet production with NLO QCD 
calculation of Un- 



expected manifestation of effects due to BFKL dynamics. Unfortunately, it 
is also the region, where playing with scales does wonders and may easily 
mask the presence of new phenomena. Although in this region the concept 
in resolved (virtual) photon does not have to be introduced it turns out to 
be very useful phenomenologically. 

From the NLO parton level Monte Carlo codes appropriate for this kine- 
matic region three are most frequently used in phenomenological applica- 
tions: DISENT m, NLO JET and JETVIP 0. Prom them only JETVIP 
includes the contribution of resolved virtual photon and only NLO JET offers 
the option of calculating three jet production at the NLO. 

Out of many interesting analyses of jet production I will briefly mention 
two. The series of investigations of dijet production in the region <C 
has shown convincingly the importance of the contributions of resolved vir- 
tual photon for the description of HERA data i7 . The relevance of the 
concept of resolved virtual photon as an approximate method for includ- 
ing higher order direct photon contributions has been studied in using 
NLOJET in the three jet mode. These studies have also demonstrated the 
importance of the contributions of the longitudinal virtual photon. 

Much attention has recently been paid to forward jet production at 
HERA in the kinematic region where Q'^/E^ is around unity. This region 
has been suggested jUj as suitable place to look for manifestations of BFKL 
dynamics. However, identification of these effects has turned out to be 
complicated by the contributions of the resolved virtual photon, as well as 
by large scale uncertainties of QCD calculations. 

The ZEUS ^Oj and HI data of forward jet production in the kine- 
matic region 0.5 < E^/Q'^ < 2, Et > 5 GeV and Pz/Ep > 0,05 have been 
analysed in NLO QCD using JETVIP and including the resolved virtual 
photon contribution jl2j . The common scale M was set to M = ^Jq'^^e\. 
and allowed to vary between M = ^/q'^+e^^ and M = ^/q^+se^. 
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Fig. 2. The HI data on forward 7r° production ^1] compared to calculations of jT^] 
(left) and [[5] (right), which include only direct photon contributions. The band 
between two dashed curves in the left part and the shaded bands at right correspond 
to the variation of the common scale in the interval •^(Q^+£;^)/4 < M < y^Q'^+E^. 



The results, presented in Fig. ^ show that direct photon contribution 
alone undershoots data significantly, but the inclusion of NLO resolved pho- 
ton contribution leads to nice agreement! However, as shown by the dotted 
and dashed curves, this agreement relies crucially on the chosen scale. 

4. Inclusive particle production in 7*p collisions 

Similarly as for jets the interest in this process is motivated primarily by 
search for BFKL effects. In jl3j the recent HI data on large px forward vr^ 
production jl4j in the region of moderate were compared to QCD calcu- 
lations including only the direct photon contribution. Setting the common 
scale fi = M = ^/(q^+e^)/2 the authors found a nice description of the HI 
data for broad range of values of (see left part of Fig. I^J. However, they 
noted strong dependence of their results on M. Varying their preferred 
value of M by factor of mere \/2 their results change by a factor of more 
than 2! This significant scale dependence in their view "suggests the pres- 
ence of non-negligible NNLO effects" and together with large X-factors " 
restraining, for the moment, any empirical suggestion of dynamics different 
to plain DGLAP evolution.". I share their reservations. 

The same HI data on forward vr'^ have been analyzed also in using 
the same input and identical choice of scales as in J^l- Not surprisingly 
they have also found a good agreement with HI data (see right plot of 
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Fig. 3. The HI data on forward tt" as a function of xbj for pT > 2.5 GeV (left) 
and of pt for three different intervals in compared to calculations of ^Hl with 
(solid curves) and without (dashed curves) the resolved photon contribution. 



Fig. 12) and factor of 2 difference between their results corresponding to 
M = ^/{Q^+El)/4 and M = ^/q^+e^. Compared to I :! the authors of [1^1 
are more optimistic and claim that their "default predictions, endowed with 
theoretical uncertainties estimated by moderate unphysical scale variations 
led to a satisfactory description of the HERA data in the preponderant part 
of the accessed phase space." In view of the arbitrariness of the standard 
choice of the common scale I consider their conclusion difficult to justify. 

This is illustrated by the most recent analysis of the same HI data per- 
formed in taking into account the contribution of the resolved virtual 
photon. This paper includes the most detailed analysis of the dependence 
on different scales, finding very different dependence on the renormaliza- 
tion, factorization and fragmentation scales. In particular "large instability 
is observed when varying independently the renormalization and fragmen- 
tation scales." Default calculations (see Fig. ^ have been performed for 
the common scale set to M = ^/q'^+e^. Compared to [2] ^-nd ^3] their 
common scale is bigger by a factor ^/2 and thus their NLO direct contribu- 
tion substantially smaller with the resolved photon contribution filling the 
gap! However, again large scale sensitivity "prevents a really quantitative 
prediction for the single pion inclusive distribution in the forward region." 

5. Jets and inclusive particle production in 77 collisions 

Inclusive particle production in 77 collisions has been a challenge for per- 
turbative QCD already at PETRA and has become so even more at LEP2. 
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Fig. 4. Left: the comparison of L3 data on inclusive charged particle production in 
77 collisions J7] with the NLO parton level calculation of JS|. Middle: the same 
data compared to the predictions of PYTHIA and PHOJET LO event generators. 
Right: L3 data on jet production compared to NLO QCD calculations of j2()| . 

The recent L3 data on charged particle production in 77 collisions jl7j . 
shown in Fig. 0J are far above the NLO QCD predictions ^S] in the range 
of transverse momenta, where the direct-direct contribution dominates and 
where there is thus little chance that higher order QCD corrections could 
make any difference. There is, however, one aspect of the L3 analysis which 
deserves closer attention. As shown in the middle plot of Fig. |1J precisely in 
the region where data exceed by more than order of magnitude NLO QCD 
prediction there is a huge difference between two LO MC event generators, 
PYTHIA and PHOJET, the former even exceeding the data! Both these 
MC use the same LO cross sections, which dominate the large pT region, 
but nevertheless lead to vastly different shape of the distribution. This rep- 
resents a warning that the effects beyond purely perturbative stage of the 
process and where the MC do differ, may play a crucial role for the tail of 
distributions like that in Fig. ^ 

The disagreement between L3 data and QCD calculations extends, as 
shown in the right plot of Fig. |1J also to jet production. Note that for jets 
both the shape and magnitude of L3 jet data and NLO QCD predictions 
of j2L)j are significantly different. However, also in this case PYTHIA lies 
above and PHOJET below the L3 data. The puzzle is still there. 

6. Summary and Conclusions 

Jet and inclusive hadron production in 7p, 7*p and 77 collisions has pro- 
vided wealth of new data for the comparison with higher order QCD cal- 
culations. Unfortunately, in the most interesting processes, like forward 
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jet and hadron production at HERA, where signals of new physics are ex- 
pected, theoretical predictions are burdened with large uncertainties stem- 
ming primarily from strong dependence on renormalization, factorization 
and fragmentation scales. The standard procedure of identifying all these 
scales and setting them equal to some physical scale characteristic for the 
process lacks theoretical justification and, moreover, is ambiguous as it im- 
plicitly assumes working in ad hoc chosen MS RS. Systematic investigation 
of scale dependence of QCD calculations should be performed with the aim 
of formulating better justified procedure for choosing the various scales. 

In theoretically seemingly clean case of jet and inclusive hadron produc- 
tion in 77 collisions perturbative QCD faces serious challenge to explain 
recent L3 data. 
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